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Abstract 

We report on the possibilities of using the method of normal funda- 
mental systems for solving some problems of oscillation theory. Large 
elastic dynamical systems with continuous and discrete parameters are 
considered, which have many different engineering applications. Intensive 
oscillations in such systems are possible, but not desirable. Therefore, it is 
very important to obtain conditions for which oscillations take or not-take 
place. Mathematically, one needs to search for the solutions of partial dif- 
ferential equations satisfying both boundary and conjugation conditions. 
In this paper we overview the methodology of normal fundamental sys- 
tems for the study of such oscillation problems, which provide an efficient 
and reliable computational method. The obtained results permit to an- 
alyze the influence of different system parameters on oscillations as well 
as to compute the optimal feedback parameters for the active vibration 
control of the systems. 
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1 Introduction 



Boundary problems with continuous and discrete parameters constitute a wide 
class of oscillation problems with a great relevance in mathematics and applica- 
tions. They appear, for instance, in the study of mechanical oscillations of ship 
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hulls, ship masts, antennas, supporting surfaces of aircrafts, turbine fins and 
shafts, structural components of automatic control apparatus, etc. The struc- 
tural study of the dynamics for such elastic systems imply solving some partial 
differential equations with non-stationary continuous and discrete coefficients. 
Moreover, the solutions of the partial differential equations have to satisfy not 
only boundary conditions, but also some conjugation constraints, which con- 
siderably complicate the problem. We claim that the method of the normal 
fundamental system of solutions is a very efficient and reliable way for numeri- 
cally solving such boundary problems with continuous and discrete parameters. 



2 The Method 

Many important problems of oscillation theory may be reduced to the study of 
systems of partial differential equations with variable coefficients, which describe 
oscillation processes in very different systems and which may be written in the 
form: 

ft N N & N f) 

^ = £ A kj{y>j + E + E , 

y 3=1 3=1 3=1 V > 

Vo < V < Vn , k = l,N, 

where Zj(y,t) is the deviation function of the oscillation system from an equi- 
librium position, and the coefficients Akj(y), Bkj(y) and Ckj{y) are real- valued 
piecewise continuous functions bounded on the interval [yo,y n ], i.e., \Akj(y)\ < 
M, \B k j(y)\ < M, \Ckj(y)\ < M for some constant M and k,j = 1, N. 
Besides, we have the next linear homogeneous boundary conditions: 



z r (0) = 0, r = l,m, 
z s (l) = 0, s = m+l,N, 

where from here always TV = 2m. 

We assume that the conjugation conditions take place at the given arbitrary 
discontinuity points y = y» of functions Akj(y), Bkj(y) and Ckj(y), and they 
may be presented as follows: 



k=i fe=i 



dt 2 



dt 

k=i fc=i 



fc=i m fc=i m 



where p — 1,N, i — l,n— 1 and the superscript i means the value of the 
corresponding function on the i-th interval. In many practical applications 
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conditions J3J are not taken into account because of significant complication of 
investigations and this causes qualitative and quantitative errors in the results 

El 

The continuous-discrete boundary problem serves as a model for a 

large number of oscillation dynamical systems 

The standard way for integration of system consists in applying the 

Fourier method of separation of variables [3| . The solution then is searched in 
the form of a product of two functions, where the first one depends only on the 
coordinate y, and the second one depends only on time t: 



z k {y,t) = ip k {y)e xt . 

The problem is thus reduced to the one of finding the values for the parameter 
A, at which there exist non-trivial solutions <pk(y) of the system Such 
values A are called eigenvalues of the boundary problem, and the corresponding 
solutions (fik(y) are called eigenfunctions. The eigenvalues for the oscillation 
system, described by equations are complex quantities A = q + ip. The 

oscillations have increasing or decreasing amplitude, depending on the sign of 
the real part q of the eigenvalue. 

In many practical oscillation theory problems, one is interested to know 
eigenvalues at which oscillations with a constant amplitude occur. That is only 
possible at an imaginary eigenvalue A = ip, where p is the eigenfrequency of 
the oscillations of the considered system. We assume that the system makes 
harmonic oscillations with some constant amplitude, searching for harmonic 
functions z^ (y, t) of time. In this case we have ifk (y) = (p^' (y) + i^p^ (y) , where 

fk iv) an< ^ fkiv) are real- valued functions. Hence, the next form of solution 
may be used: 

Zk(v,t) = [<p£\v)+i$\v)]e* t . (4) 
In particular cases, one of the functions ip^ (y) or ip^jp (y) may be identi- 

(2) 

cally equal to zero. Thus, if <p k (y) = 0, then fk{y) is a real- valued function. 
This happens, for instance, in oscillatory motions of mechanical systems with 
no forces proportional to the velocity, i.e. in the absence of the terms with 
d z 

derivatives — 
dt 

After substitution of into equations Q and separation of imaginary and 
real parts, we obtain respectively: 



, (l) N N N 



dy 



J = l 3=1 3=1 



, (2) N N N 

=E A «(f)^ 2) (f) -p 2 Y, B v(yM 2) (y) . ( 5 ) 

aV 3 = 1 3 = 1 3 = 1 



where fc = 1, N. 
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The boundary conditions @ after similar transformations may be written 

as 

^\y )=<p?\y Q )=Q, r = l~^, 
<P?Hv n ) = <P?Hv n )=0, s = m + l,N. 

Analogously, the conjugation conditions J2J will be transformed to the next 
ones: 



N N 

fc=l fc=l 

AT iV 

- p E d £J 1} (») = E a k l ] P (vi)Vkl & 

k=l k=l 
N N 



- p 2 E O^S (w) - pE 4%)*$ (»*) > 
fe=i fc=i 

iV iV 

E (») - p 2 E (») 

fc=l fc=l 

AT iV 

fe=l fc=l 

N N 

-p 2 E O^S m+pE d£(w)*S (W) > (7) 



fc=l k=l 



where p = 1,N. 

Since the initial problem is reduced to the form ©-0, we can con- 

sider now the boundary problem for the system of ordinary differential equa- 
tions, or, more precisely, the problem with the normal system of linear ordinary 
differential equations: 

, N 

d?=I>y(w)*i> k = T^, (8) 

V 3=1 

where a,kj(y) are piecewise-continuous functions on [yo,y n ] satisfying Lipschitz- 
type conditions and, therefore, bounded, i.e. \dkj{y)\ < M, with constant M 
and k,j = 1, AT. 

We have for system |JSJ the next N linear homogeneous boundary conditions: 



z r {yo)=0, r = 1,771, ^ 
z s (y n ) = 0, s = m+l,N. 



At the discontinuity points y = y^ i = 1, n — 1, of functions akj(y), the 
solution Zfe(y) of the system (JHJ should satisfy the next conjugation conditions: 
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N 



N 



XX.fe (?/*) 4 (Vi) = £ b j,k iVi) Z k iVi) > 



(10) 



fe=l 



j = 1, N, i = 1, ?i — 1, and where d 3 ^{yi) is the value of the function at the ?-th 
interval, and bj^{yi) is the value of the function at the (i + l)-th interval. 

Now, we consider a Cauchy problem for the system © for each interval 
Vi-i < V < Hi, i — 1, n, with the initial conditions at points y — y%—\. According 
to Picard's Theorem, each such problem has a unique solution satisfying the 
given initial conditions. Therefore, we are in conditions to apply some well- 
known numerical method, for example the Runge-Kutta method, and find the 
fundamental system of solutions for each i-th interval, i.e., to find N linearly 
independent solutions of the system (JHJ) in the form 



$W (y) 



= Wk,j{y)}k,i 



(ii) 



where each function (f^j(y) is defined and continuous at i = l,n; k 

is a solution number and j is a function number, k,j — l,N. At the points 
y = yi-i, these solutions satisfy the conditions 

^(y i . 1 ) = E, i = T^, 

where E is the unit matrix, or, in a more detailed way, 




i = 1, n. 



(12) 



Such fundamental system of solutions (|llf> with initial conditions i|12|) is called 
a normal fundamental system 0. We remark that for each point y — yi-\, 
i = 1, n, the existence and uniqueness theorem applies. The general solution of 
the system (JHJl at the i-th interval may be written with the help of the normal 
fundamental system of solutions as 



N 



l%) = £q i V2(y), 

4=1 



(13) 



with Cj arbitrary constants and yi—i < y < y%, k — 1, N, i = 1, n. Satisfaction 

of the conjugation conditions (|10|l allows to obtain the coefficients in the 
form 



N 



c (i) = V C {1) u (l) 

g=m+l 



(14) 
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with the coefficients v/p q given by recurrence as follows: 



,(*+!) _ 



N 



N 



N 



fc=l 



iV AT ... ... 

b N,j-i E rf jv,fc E u p,Wpi(2/i) b Nj+ i 

k=l p=l 



b N 



N 



bN,N 



^ (15) 

where j = 1, N, q = m + 1, JV and i = 1, n — 1. 

For application of the formulas Ijf 5|l we need to know the initial coefficients 
il for all j = 1,N and q 
V = Ho we can find 



m+ 1, JV. From the boundary conditions JjJJl at 
1, 3 = q, 



"SHn ' *' j = l,N, q = m + l 1 N 
0, j^q, 



(16) 



Thus, the solution of the system © for y Q < y < y n , which satisfies to the 
conjugation conditions (JTDJ , is determined by the formula i|T3|) subject to l{Tl)l- 
l|f 61) . According to the existence theorem, such solution exists and it is unique 
at the interval yo < y < y n . 

Imposing also the initial condition tftfy at u = tu,. we obtain a homoge- 
neous system of m linear algebraic equations in the coefficients C,p for all 
q = m + l,N: 

N N 



where s — m + I , N. 

The boundary problem (f5|) — (| 1 01) has a non-trivial solution only at the con- 
dition of existence of non-trivial solution of the system (|f 7|) . 

Thus, we conclude that a necessary and sufficient condition for the existence 
of a non-trivial solution of the boundary problem JHJ)-(|!5J) is given by 



D = det 



JV 



= 0, 



(18) 



s, q = m + 1, AT. Taking into account the dependencies of the functions f^liyn) 
on the value A, we determine the eigenfrequencies of the oscillation system as 
roots of the equation (|I8|I . 

Given a concrete system, the necessary and sufficient condition l|I8|) provides 
a general method to investigate the frequency spectrum of possible oscillations 
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and its dependence on different parameters of the system. In the next section 
we review the application of this method for the study of different oscillation 
problems. 



3 Some Applications 

Let us mention in this section some applications of the method of normal funda- 
mental systems of solutions for solving oscillation problems in different technical 
systems. 

The authors of apply this method for analysis of dynamics of some con- 
trolled machine units, which are widespread in different modern apparatuses, 
machines for materials processing, transportation, etc. The considered con- 
trolled machine units consist of the motor, some mechanical element and the 
elastic connection of significant length between them. The mathematical model 
of such controlled dynamical system with distributed and discrete parameters 
is a partial differential equation satisfying non-classical boundary conditions as 
follows 

„ d 2 ip d 2 ip d 3 ip , , 

GI ^-^W + C ^ = °' (19) 

d 2 ip d 2 ip dip „ dip , . 

d 2 ip d 2 ip dip dip 

* = L > j2 ^ + ^di^ + GI ^~ m ^m=°> (21) 

where ip(x,t) is the shaft arbitrary section twist; G is the shear modulus; I p 
is the polar inertia moment; p is the line density of the unit of the elastic 
link material length; J\ is the inertia moment of the motor rotor; J2 is the 
reduced shaft inertia moment of intermediate mechanisms of the motor and the 
tool; £1 is the dissipation coefficient in the shaft; /3 is the rigidity of the motor 
mechanical characteristic; mo is the tool mass; ax is the coefficient of the linear 
rate component of the tool resistance moment. 

The numerical method of the normal fundamental systems allows to deter- 
mine the frequencies of possible oscillations for the case of the problem (flfjp — . 
The influence of different parameters on the frequencies is also analyzed by the 
authors. It is found the possibility to avoid undesirable oscillations with the help 
of application of different feedbacks, which purposefully change the rigidity of 
the motor mechanical characteristic in the machine unit. 

In the paper 0, we apply the numerical method of normal fundamental 
systems of solutions for investigation of oscillations in a flexible clastic system 
of a spacecraft of big size. 

The basic mathematical model, which is used for investigations, is: 
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x = : u (0, t) = , 



(23) 



x = l: mESP——+m{d + b)— + ES(m + l3b) 



dxdt 3 dt 3 dxdt 2 

where w(x,f) is the displacement of a point of the bar with abscissa x, S the 
bar cross-section area, p the unit volume mass, E the modulus of the material 
elasticity, (3 the coefficient of dissipation in the bar material, b the damping 
factor of the executive mechanism, c the rigidity of the centering spring, d the 
feedback coefficient, m the specimen mass, and I is the length of the elastic bar. 

The authors show that for each such dynamical system it is possible to char- 
acterize possible oscillations using normal fundamental functions. Moreover, 
it is possible to change, purposefully, eigenvalues of the considered nonconser- 
vative boundary problem. With this aim, different feedbacks may be applied, 
changing the feedback coefficient d in a desirable way. This allows to avoid in 
the system an excitation of undesirable oscillations, including self-excited ones. 

These results of the considered numerical method coincide with the ones 
found in obtained to the problem 122|) - l|24|l with a different approach: the 
approximate analytical method of small parameter. 

The analysis of dynamics for some aerial cable- way systems, by the method 
of normal fundamental systems of solutions, may be found in 0. If discrete 
loads with masses mi, TO2, m„_i move along the cable- way at a steady 
speed v, then the mathematical model for such controlled dynamical hybrid 
system is the next non-stationary boundary problem: 

d 2 Ui d 2 Ui 

P ^F = T ^' (25) 

Ul (0,t) = u n (M) = 0, (26) 

with conjugation conditions 

u i+1 (k, t) = Ui(k, t)\ x=h ( t) , (27) 



,d 2 u l d 2 Ul 2 52M m .p/d^+i dui 

where t is the time variable; x is the coordinate of some point of the cable; p 
is the mass per unit length; T is the tension in the cable; I is the length of the 
cable; mi are the masses of the discrete loads at arbitrary points x = U (t) for 
all i = l;n— 1; li(t) = U + vt, lo(t) =0, to < t < t\\ h is the coordinate 
of the arbitrary point with the discrete mass rm at the initial moment of time 
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t = 0, i = 1, n — 1, Iq = 0, l n = I; Ui(x, t) are the functions of deviation from 
equilibrium position for the point x € at the moment t, i = l,n. 

For the problem (|25|) - (|28|l . which is a continuous-discrete problem, the 
method of the normal fundamental functions is very efficient. It allows to find 
the dependence of the dynamical characteristics of the system on the controlled 
parameter, which is the speed v. 

Finally, let us mention one open problem which may be solved by the method 
of normal fundamental functions. In the paper g], the approximate analytical 
method is used for investigations of possible oscillations in the transmission 
pipelines intended for lifting minerals from great depth. 

A pipe of significant length with a big mass platform on the end is a basic 
constructive element of such deep-water plants. Intensive dynamical processes 
of different physical nature may be excited under the influence of forced wave 
disturbances and nonlinear hydrodynamical forces in the system. Water flow, 
which washes the pipe, can continuously transport the energy into the system, 
significantly exciting the system and causing dangerous oscillations. As a result, 
breaking of the plant is possible. 

The mathematical model of vibrations excitation, in the case when the pipe is 
washed by a water flow, is a dissipative wave equation of longitudinal vibrations: 



2 u „ d 2 ( „du 



a 2 



= 0, (29) 



dt 2 dx 2 \ ' dt, 

E 

where u(x,t) is the longitudinal displacement of pipe points; a 2 = — ; E is the 

P 

modulus of elasticity of the material; p is the mass of material volume unit; j3 
is the coefficient characterizing internal friction in the material. 

For transmission pipelines with the clastically attached pipe, boundary con- 
ditions for the equation (|29|) are the following: 

x = 0, ES^(u(0,t)+p^^\=ku(0,t), (30) 
#HM) LT?q d ( a du(L,t) 

du(L, t) f du(L,t) \ 3 
= ai ^^- a3 {^^) ■ (31) 

where M is the platform mass; S is the pipe cross-section area; k is the lon- 
gitudinal rigidity of the elastic hanger at the point of pipe attaching; L is the 

. , , du(L,t) fdu(L,t)\ 3 . , , 
pipe length; at\ a 3 is the force of the non-hnear resis- 

ot \ at J 

tance, which appears as a result of the interaction between the body and the 
environment; a%, «3 are constants characterizing the resistant environment. 

It is obvious that the method of normal fundamental functions may be used 
for numerical solving the boundary problem l|29|) - (|31|l . which was solved in the 
paper 0] only approximately and with some additional assumptions. 
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For all the above mentioned systems, the technique for determining complex 
eigenvalues of the corresponding boundary problem is developed. With the help 
of the method of normal fundamental systems of solutions, the dependencies of 
the oscillation frequencies on different physical parameters of the systems are 
obtained. Some possibilities of control of the frequency spectrum, in which 
oscillations are possible, are found for each mentioned system. 
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